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We study a thermally induced spin flip of an electron spin located in a semiconductor quantum 
dot. This interesting effect arises from an intriguing interplay between the Zeeman coupling to 
an external magnetic field and the hyperfine interaction with the surrounding nuclear spins. By 
considering a minimal model, we explain the main mechanism driving this spin flip and analyze 
its dependence on the strength of the external magnetic field, the number of nuclear spins and the 
ratio of the electron and nuclear Zeeman energies, respectively. Finally we show, that this minimal 
model can be applied to experimentally relevant QDs in III-V heterostructures, where we explicitly 
predict the temperature at which the spin flip occurs. 
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I. INTRODUCTION 

The physics of an electron spin confined to a quantum 
dot (QD) has been in the focus of condensed matter re¬ 
search for many years due to its possible applications in 
quantum information processing. Following the seminal 
work 1 of Loss and DiVincenzo, the two eigenstates of the 
electron spin can be used to define a quantum bit (qubit). 
However, in the host nanostructures, the e lect ron spin is 
subject to a variety of different interaction^^ leading to 
a loss of the stored information and decoherence. Im¬ 
portant examples for these interactions are the hyperfine 
interaction (HI) with surrounding nuclear spins or the 
coupling of the electron spin to electrical fields via spin 
orbit interaction. The dynamics of the electron spin in¬ 
fluenced by these effects - alone or combined - have been 
studied in great detaiP®31 j n the last years enhancing 
our knowledge about spin physics in QDpHm. Along 
with these insights also strategies were developed to re¬ 
duce the loss of information-®. However, while electrical 
noise seems fairly controllable, the HI with the nuclear 
spins remains nonetheless as a major source of decoher¬ 
ence, since nuclear spins are intrinsic to most QD host 
materials such as the widely used III-V heterostructures. 
Yet recently, also a change of paradigm can be observed, 
where the unavoidable interaction with the nuclear spins 
is considered rather as a resource for interesting physics 
than an obstacle!® 

In this notion, we want to present our findings on a 
thermally induced flip of the electron spin. These results 
were obtained by applying standard statistical physics to 
a minimal model for spin dynamics in a QD consisting 
of an external magnetic field to which the electron and 
the nuclear spins couple and the HI which links the elec¬ 
tron spin to all nuclear spins. The effect of an external 
magnetic field on the electron is described by the Zeeman 
Hamiltonian 

Me = 9*PbB z S z , ( 1 ) 

where g* is the effective g-factor 20 of the electron, is 
the Bohr magneton, and S z is the electron spin compo¬ 


nent parallel to the external magnetic held B z . Likewise, 
nuclear spins also align with respect to the external mag¬ 
netic held by means of another Zeeman term 

K K 

Bze = X = ~ X 9NPNB z ik,z j (2) 

k= 1 fc=l 

where gw is the nuclear g-factor, /rjy is the nuclear mag¬ 
neton and the sum is running over all nuclear spin z com¬ 
ponents Ik , z • Note the relative sign difference between 
Eq. ([I]) and Eq. (|2j) , which resembles the negative sign of 
the electron charge. 

By writing the nuclear Zeeman Hamiltonian in the 
form of Eq. ([2]), we already assume that there is only 
one spin species present in the dot. This assumption 
will simplify our reasoning and, thus, allows us to better 
identify the main physics being relevant for our results. 
Distinguishing different spin species would not change 
our hndings signihcantly as we discuss in the end of this 
article. This specihc choice also simplifies the HI Hamil¬ 
tonian 


K 
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H m = A m X 1‘KuOI 2 [Szik,z + §( S + I k + S_I*,+) 

k =1 


(3) 

where \4>(rk)\ 2 is the probability to find the electron at 
the site of the fc-th spin-carrying nucleus. The energy of 
the HI is given by a constant 


Am = 9 ■ 9n ■ C 


( 4 ) 


with C > 0 being a material dependent energy scale. 
Since the HI interaction is strongly localized around the 
res pectiv e nucleus, the bare electron g-factor g sa 2 enters 
here® 21 . Thus, the sign of the HI is determined by the 
sign of the nuclear g-factor gw- This sign, however, plays 
an important role, since it determines the form of the 
ground state of the HI. 

If the coupling constant Ahi is positive (negative), the 
ground state of the bare HI will favor an anti-parallel 
(parallel) alignment of the electron spin with respect to 
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all nuclear spins. Both ground states are twice degener¬ 
ate, since a flip of all spins results in the same energy. 
Similarly, the two Zeeman terms also show the same two 
types of spin ordering depending on the signs of the re¬ 
spective g-factors g* and g^. In contrast to the HI, these 
ground states are unique. If, for instance, g* < 0 and 
gN > 0, the Zeeman terms would force both the elec¬ 
tron spin and all nuclear spins to be parallel to the ex¬ 
ternal magnetic field at zero temperature. Thus, when 
both the HI and the Zeeman interaction are present, 
there can arise an interesting competition of spin order¬ 
ing with the electron spin being parallel or anti-parallel 
with respect to the nuclear spins. In particular, if the 
external magnetic field is sufficiently small, the HI is still 
strong enough to maintain the anti-parallel alignment of 
the electron spin with respect to the nuclear spins. If 
then additionally the signs of the g-factors are given by 
g* < 0 and g^ > 0, the electron spin will be anti-parallel 
to the external magnetic field, whenever its Zeeman en¬ 
ergy is below the total Zeeman energy of all nuclear spins. 
Starting from this particular ground state, we will show 
that a sudden flip of all spins can happen at a finite tem¬ 
perature To > 0. 

The article is organized as follows: In Sec. II.A, we 
will argue how the HI Hamiltonian can be simplified 
based on physical arguments. In Sec. II.B, we will then 
apply standard statistical physics to the Hamiltonian 
H = H% E + H% e + H m , where we neglect the off-diagonal 
parts of the HI. Doing so, we calculate the thermal ex¬ 
pectation value of the electron spin, whose properties are 
studied both analytically and numerically in Secs. II.C 
and II.D, respectively. After this mathematical analy¬ 
sis, we will then explain the physical mechanism being 
responsible for this spin flip in Sec. II.E. In Sec. Ill, 
we finally review our initial simplifications of the HI and 
discuss in which real systems our findings should be ob¬ 
servable. In Apps. [A] and [Bj we calculate the thermal 
expectation values for the full HI Hamiltonian including 
its off-diagonal part, i.e. the flip-flop terms. In App. 
[Cl we analyze the behavior of the electron spin at zero 
temperature for this interaction. 

II. THERMAL ELECTRON SPIN FLIP 
A. Simplified Hamiltonian 

As mentioned above, we will first introduce certain 
simplifications to the HI Hamiltonian, which allow for 
an analytical calculation of thermal expectation values: 

1. We assume that all nuclear spins are spin one-half, 
where the number of nuclear spins is K. 

2. We will use the so-called box-model, where the 
probabilities |</>(r/)| 2 = 1 /K are all the same. By 
this, we assume that every nucleus in the dot car¬ 
ries a spin and that the envelope function <j>{fk) of 
the electron does not change much inside the QD. 


With these two assumptions and g* < 0 as explained in 
the introduction, the total Hamiltonian is given by 

H a = e + -Hze + Hhi = 

- \g*\g,BB z S z - gNHNB z J z 

+ -^-[S Z J Z + — (S+J_ + S-J +)], (5) 

where we have introduced the total nuclear spin J = 
Y2k=i Ik f° r convenience. The usual raising and lowering 
operators S± = S x ± i S y and J± = I x ± i J y form the flip- 
flop terms. By means of the parameter a we distinguish 
between the full Hamiltonian (a = 1) and a simplified 
Hamiltonian (a = 0), which allows us to present the basic 
physics of the electron spin flip more easily. 

Before we proceed with the calculation of thermal ex¬ 
pectation values, we will introduce dimensionless units by 
measuring all energies in units of The total Hamil¬ 
tonian then reads 

H a = —<jS z - vJ z + 2S z Ik, z + a{S + J- + S-J +), (6) 

where the dimensionless parameters are given by a = 
I 9*\hbB x /4$ and v = g N PNB z /j^. These two 
quantities are obviously not independent of each other 
since both are proportional to the external magnetic 
field B z . Thus we choose a = Kpv , where p = 
\g*\p,BB z /(Kg^ppfB z ) is the ratio of the Zeeman en¬ 
ergies of the electron and all nuclear spins, respectively. 
This ratio can be also characterized by the critical num¬ 
ber k = pK , which is a constant for a given material. 

In the following, we will first analyze the thermal ex¬ 
pectation value of the electron spin for the simplified 
Hamiltonian Hq. Since this Hamiltonian is already di¬ 
agonal in the basis of products states of the individual 
spin states, all calculations are much easier and, hence, 
the physics causing the spin flip becomes more apparent. 
However, it is not clear in the first place, if the neglected 
flip-flop terms give rise to quantum fluctuations, which 
destroy the electron spin flip. Thus, we also have calcu¬ 
lated the thermal expectation value of the electron spin 
for the full Hamiltonian Hi including the flip-flop terms. 
The details of these calculations are reported in Apps. 
0to0 while we will only present the respective results 
in the main text. Interestingly, many findings are un¬ 
changed with respect to the simpler case or restored in 
the limit of large system sizes, where the flip-flop terms 
are shown to be irrelevant. 


B. Thermal expectation values 

Without the flip-flop terms, the simplified Hamilto¬ 
nian is already diagonal in the product basis |ms) <8> 
<S>k=i \ m k), where | ms) is an eigenstate of S z with 
ms € {—1/2,1/2}. Similarly, the state \mk) is an 
eigenstate of the fc-th nuclear spin operator /fc i2 with 
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m k £ {—1/2,1/2}. Thus, the partition function for this 
Hamiltonian is easily calculated in the product basis 


Z = Ti[e-^ o/kBT } = 

E E 

m S {rn k }£ =1 


E E e -(ms,m K ,m K -i,...\H 0 \m s ,m K ,m K -l ,■■■)/k B T 


(7) 


where the diagonal matrix elements are given by 


- {ms,m K ,m K -u ... |flo|m s , m K , m K -i ,.. .)/k B T 
„ K 


= - [pisK ms + (ism k - 2 msm k )\ 


( 8 ) 


k =1 


with the dimensionless temperature r = fcsT/^gt. Ex¬ 
ploiting the fact that sums in the exponential functions 
factorize finally yields 


K 


z - e n e en 

ms k =1 mfc 


= (ei i + ei _i) A + (e_i i + e_i _i) K , (9) 

' 2 ’ 2 2 ’ 2 7 v 2 ’ 2 2 ’ 2 7 v/ 


where 


e ms ,m fc = exp[— (pism s + ism k - 2m s m k )} ■ (10) 

T 

With the partition function at hand, the calculation of 
the thermal expectation value of the electron spin is read¬ 
ily obtained 

« t d t d 

(S Z ) T = -^ln [Z] = — i-ln[Z] 

V 2 da 2Kis dp [ 1 
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where the function 

n = n(r, v, p, K) = 


e_ i i + e_ i _ i 

2 ’ 2 2 ’ 2 


e i i + ei _ i 
2 ’ 2 2 ’ 2 


K 


(ii) 


( 12 ) 
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FIG. 1: (Color online) a): Yl(r,is, p, K) as a function of tem¬ 
perature r for v = 2, p = 0.35 and A' = 10 4 . As we show 
in the text and in Fig. [2j this choice of parameters fulfills the 
necessary conditions for the spin flip. At a temperature to, 
the function suddenly drops from very large values to zero. At 
temperatures r > n the function surprisingly rises again and 
saturates at II = 1 for large temperatures, b): The thermal 
expectation value of the electron spin ( S Z ) T exhibits a sud¬ 
den flip at to. For temperatures above ti the electron spin 
is thermally equilibrated, c): The thermal expectation value 
of the total nuclear spin ( J Z ) T = j/y In \Z\ decreases before 
the electron spin flip. At to it is suddenly increased again. 
For v < 1 instead, the total nuclear spin would exhibit a flip 
similar to the electron spin. 


C. Analytical analysis of the spin flip 

By inserting the definition of the exponential functions 
in Eq. (10) into Eq. (12) and rearranging factors, we find 


ir K = U(t,is iPi K) = 


controls the behavior of ( S Z ) T . If n —>• oo, we find 
(S Z ) T —1- —1/2, n = 0 results in (S Z ) T = 1/2 and, fi¬ 
nally, n = 1 yields ( S Z ) T = 0. As it turns out, the 
electron spin has to go through exactly these steps for 
the thermal spin flip to occur as illustrated in Fig. [l] In 
the following, we will first explore the parameter space of 
n(T, v , p , K) to find mathematically the conditions neces¬ 
sary for the spin flip to occur. Afterwards we calculate at 
which temperatures n and, consequently, ( S Z ) T undergo 
their characteristic changes. To this end, we will analyze 
n analytically and compare the findings with numerical 
calculations of ( S Z ) T . Finally, we will interpret these re¬ 
sults in order to understand when and why this effect 
appears in a physical system. 


exp[-] 


2 pis, exp[i{-^ - 1}] + exp [±{is + 1}] 


t exp[i{-^ + l}]+exp[i{i/-l}] 


K 


(13) 


By further rearrangements of factors in Eq. |13j) , we iden¬ 
tify that p < 1 and 0 < is < p~ x are necessary conditions 
for n to diverge at r —» 0 and, consequently, for the spin 
flip to occur. Within this parameter regime, we want to 
identify the temperature To, at which n drops from in¬ 
finity to zero, and the temperature T\ at which n rises to 
1 as indicated in Fig. |T] 

The latter temperature can be readily read off from 
Eq. (13), since n = 1 for all temperatures r well above 


Ti = 2pisK. 


(14) 
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At a specific temperature r the function 7 r defined in 


Eq. (13) changes from 7r > 1 to 7r < 1. Since one has to 


take it to the power of K 1, this marks the tempera¬ 
ture, at which the sudden drop from II » 1 to II <C 1, 
and, hence, the spin flip occurs. Thus, the transcendental 
equation defining tq reads 


r ° 1 

p = ~Tv' n 


ex P[^{~^ ~ 1}] + ex P [^{^ + 1}] 
ex P[^{-^ + !}] + exp[^{^ - 1}] 


(15) 


This equation can be expanded in powers of A <g; l ; 


-+o(4) 

T 0 T§ 


(16) 


As a consequence, the temperature tq ~ p~ l is indepen¬ 
dent of v for p <C 1. Since p = n/K is a constant for a 
given QD, this is a rather intriguing result. This constant 
being p <C 1 corresponds to a situation, where the total 
nuclear Zeeman energy is much larger than the electron 
Zeeman energy. 

Before we give a detailed physical interpretation of our 
results, let us compare these analytical results for the 
simplified Hamiltonian Hq with i) a numerical analysis 
of Eqs. © and ( |12[ ), respectively, and ii) the behavior 
of the electron spin for the full Hamiltonian Hi calculated 
in Apps. 0 to0 


D. Numerical analysis of the spin flip 


In order to verify our analytical results for Hq, we 
show density plots of (S Z ) T = (S z ) t (t,u, p, K), where 
we choose the number of nuclear spins to be K = 10 4 . 
With K being fixed, r, p, and v remain as parameters. 
In Fig. [ 2 ] a), we show ( S Z ) T as a function of r and p 
with v = 2. As we will show below, for this choice of 
v, the electron Zeeman energy competes with the total 
HI energy. If pv > 1 the Zeeman energy exceeds the 
HI energy and the electron spin is up for all r, which 
is clearly shown in Fig. 2] a). Additionally, we plotted 
Eqs. (14) and (15) in order to demonstrate the behavior 


of 7~i and To, respectively. Both analytical results show 
a remarkable agreement with the numerical findings. In 
Fig.|b), we show the same plot for v = 0.09. For this 
choice of v, the electron Zeeman energy is always smaller 
than the total HI energy. However, if p > 1, the electron 
Zeeman energy is larger than the total nuclear Zeeman 
energy and, consequently, the electron spin flip is up for 
all r as can be seen from Fig. [ 2 ] b). 

As indicated above, in a real system p is rather a fixed 
parameter than a real variable. Hence, we also calculated 
( S Z ) T as a function of temperature r and the Zeeman 
energy of a single nucleus v, which is proportional to 
the magnetic field B z . The result is shown in Fig. [ 2 ] c), 
where we chose p = 0.09. In this figure, the behavior of 
r 0 = p _1 is most prominent. Moreover, it is obvious that 
pu < 1 is indeed a necessary condition for the spin flip. 


Finally, we have confirm that the influence of the flip- 
flop terms, which are present in the full Hamiltonian Hi 
does not destroy the electron spin flip. As we show in 
Apps. |A|an m the thermal expectation value of the elec¬ 
tron spin can be exactly solved for the box-modelP^HISI 
By investigating the temperature dependence of ( S Z ) T for 
up to K = 60 nuclear spins, we see that both tempera¬ 
ture scales To and t± are unchanged. For small numbers 
of nuclear spins, however, we find that the additional 
interaction alters the behavior of (S Z ) T . To be more spe¬ 
cific, the minimum of the thermal expectation value (S z )0 
at zero temperature is larger than —1/2 for few nuclear 
spins as can be seen from Fig. [4j For larger system 
sizes, the original value of {S z )q = —1/2 seems to be 
restored. Identifying the ground state of the system in 
App. § we are indeed able to show that (S z )0 = —1/2 
is exactly reached in the limit of large K. Moreover, the 
maximum of the electron expectation value at approxi¬ 
mately t = 2/p is altered by the flip-flop terms as is ob¬ 
vious from Fig. [5] Again, their effect is most pronounced 
for small K, while the results for Hq are reproduced for 
large system sizes. Thus, the quantum fluctuations do 
not destroy the spin flip. Physically, it seems very likely, 
that the nuclear Zeeman energy additionally stabilizes 
the spin system against the flip-flop terms. In the limit 
of large system sizes, the physics of the simplified Hamil¬ 
tonian is restored, which can be understood by analyzing 
how different states are affected by the flip-flop terms. 
The states being most efficient for this interaction are 
of the form | — 1/2,1/2,..., 1/2) and vice versa, since 
the electron spin can flip with every nuclear spin. The 
states | —1/2,1/2, —1/2,1/2, —1/2,...) (and all permuta¬ 
tions of the nuclear spins) are the less affected ones. For 
large system sizes, however, the statistical weight of the 
latter states is much higher than for the former states, 
explaining the vanishing influence of the flip-flop terms 
for increasing K. 


E. Physical interpretation of the results 

So far, we have mathematically clarified for which pa¬ 
rameters one finds the spin flip. In the following, we want 
to explain why this spin flip occurs and how the condi¬ 
tions found above can be interpreted physically. There¬ 
fore, we will have a closer look on the energies of the 
Hamiltonian Hq given in Eq. Since the Hamiltonian 
is invariant under the exchange of two nuclear spins, only 
the total nuclear angular momentum Mj = Y2k is r el- 
evant resulting in N'^ {j = (k^ M j ) equal energies of size 

E™ J S = -Kpvms - Mjv + 2 m s Mj . (17) 

A plot of E™J and E^f* is presented in Fig. 3 where a 
small system size of K = 10 nuclear spins has "been cho¬ 
sen for practical reasons. Although the spin flip is not 
perfect for such a small number of nuclear spins as is obvi¬ 
ous from Fig. [ 5 ] the physics of the spin flip still becomes 
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FIG. 2: (Color online) a): The thermal expec tatio n value of the electron spin as a function of temperature r and p for v = 2. 
The green dashed line shows the defining Eq. (15 I of to. For small p one finds to = p _1 in agreement with the plot. The light 
blue line is given by ti in Eq. 


For pv > 1 the spin flip is absent. For temperatures above ti, the electron is 


( |14| ) . For pv > 1 the spin flip is absent. For temperatures above tj, the electron is thermally 
equilibrated and, hence, (S Z ) T = 0. b): The thermal expectation value of the electron spin as a function of temperature t and 
p for v = 0.09. Clearly, the spin flip is absent for p > 1. c): The thermal expectation value of the electron spin as a function 
of temperature t and v for p = 0.09. The light blue line is given by ti in Eq. (141 . The horizontal green line corresponds to 
v = p -1 and the vertical green line to t = p -1 . 
again thermally equilibrated. 


clear. The ground state is given by ms = —1/2 and 
Mj = K/ 2, whose energy eigenvalue E K y, 2 is not degen¬ 
erate. If one follows the energies E M A< 2 and E^A starting 
from Mj = K/2 , one finds that E M y 2 increases much 


faster than Ey 2 , which is also obvious from Eq. (171. 


Since the degeneracy of the corresponding energy leve 
is strongly increasing, many states with ms = 1/2 
become thermally available for finite temperatures. Once 
the temperature reaches To, there is a strong imbalance 
between the number of states with ms = —1/2 and the 
number of states with ms = 1/2, which finally causes the 
sudden spin flip. If the temperature is further increased 
above ti almost all states are reached and, hence, one 
finds ( S Z ) T = 0. 

Finally, we would like to interpret the mathemati¬ 
cal conditions on the parameters physically. Let us 
start with the constraints on p and pv. Since the p = 
\g*\pBE z l{Kg N p N B z ) is given by the ratio of the elec¬ 
tron Zeeman and the total nuclear Zeeman energy, the 
ratio being p < 1 implies that the Zeeman energy of all 
nuclear spins exceeds the electron Zeeman energy. Sim¬ 
ilarly the product pv = 2\g*\p B B z /Am tells us that the 
Zeeman energy of the electron has to be smaller than the 
total HI energy. This imposes an upper bound on the 
external magnetic field 



FIG. 3: (Color online) Energies of the Hamiltonian as a func¬ 
tion of the total nuclear angular momentum Mj = rrik 
for electron spin up (ms = 1/2, red) and down (ms = —1/2, 
blue). The other parameters are v = 2, p = 0.35, and I\ = 10. 
The width of a column is proportional to the degeneracy of 
the respective energy level given by N&j = ( k / 2 ^-m,)- The 
green and light blue line show the temperatures to and ti, 
respectively. 


temperature 


T 0 = k 


-1 

B 


Am gNfJ’N 
2 \g*\vB ' 


(19) 


B z < U max = • ( 18 ) 

2|<? I Mb 

But the magnetic field has additionally to be large 
enough in order to separate the two temperatures To and 
T\. For small magnetic fields, the former temperature 
is given by tq = p , which corresponds to an absolute 


The latter temperature ti corresponds to the Zeeman 
splitting of the electron since 

Ti =kz 1 \g*\p B B z . ( 20 ) 

For Tq < Ti the spin flip is present, which is the case if 
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9* 

9n 

K, 

Ahi 

^max 

-^min 

To 

Ti 

Material 



[10 3 ] 

[peV] 

[mT] 


[mK] 

[mK] 

GaAs 

-0.4 

1.8 

0.44 

84 

1700 

3 

2 

5.6 

49 

CdTe 

-1.8 

-6.7 

4.9 

-34 

170 

1 

2 

0.04 

20 

Gain As “ 

-4.4 

3.1 

2.6 

93 

180 

3 

2 

1.0 

54 

Tn As 

-15 

3.5 

7.9 

98 

57 

3 

2 

0.36 

57 


“Gao.47lno.53As 


TABLE I: Relevant materials and their parameters. The val¬ 
ues of g* are taken from Ref 1201 The magnetic moment g n 
and the HI constant Ahi are averaged by grjv = JT nig \7 and 
Ahi = X/ n iAm, where m is the natural abundance of iso¬ 
tope %. The values are taken from Refs. I2I3I26I The temper¬ 
atures T\ and To are calculated at B z = 0.1-B max , at which 
To = kg 1 ■ [4/(/ + l)/3] ■ is valid. By this choice, 

we also take into account that larger nuclear spin quantum 
numbers I increase To. For simplicity, we took the smallest 
quantum number 7 m j n of different isotopes present in the dot. 


the magnetic field obeys 
9nPn 


B z ^min 


iffl 


= - B 
It* Its %*Imb « 


( 21 ) 


Hence, our results still hold if the HI constant Ahi is 
replaced by ni Ahi ■ 

Besides approximations concerning the nuclear spins, 
we also simplified the physics of the electron spin by using 
the box model for its envelope function: |</>(r/)| 2 = A' -1 . 
In reality, the probability to find the electron, should 
strongly decrease with the distance from the center of 
the QD, which is often described by a Gaussian enve¬ 
lope function \(j>{fk)\ 2 oc AT -1 exp[— r k /R], Our results 
should be modified in this case by two aspects: Nuclear 
spins with \4>{fk)\ 2 A" -1 couple only very weakly to 

the electron and can, thus, be neglected. Effectively, 
this reduces the number of involved nuclear spins from 
AT to AT e ff < K. For nuclear spins in the center, one 
finds 0(\<j)(fk)\ 2 ) = AT -1 . As a consequence, the non- 
uniform HI will (slightly) lift the degeneracy of the en¬ 
ergies in Eq. and Fig. [3j but it will not change the 
energy spectrum in principle. Therefore, the main phys¬ 
ical mechanism stays the same and our results still hold. 

Finally, we have investigated the behavior of the elec¬ 
tron spin expectation value for the full Hamiltonian in¬ 
cluding the HI flip-flop terms in Apps. [A]to[C] Doing so, 
we have confirmed, that the relevant temperature scales 
are essentially the same and that the physics of the sim¬ 
plified Hamiltonian is reproduced in the limit of a large 
number of nuclear spins. 


III. DISCUSSION AND OUTLOOK 

All of the above results were obtained applying certain 
approximations, which are in general not fulfilled by real 
systems. In the following, we will discuss all approxima¬ 
tions one by one and analyze, how a more realistic model 
would change our results. 

First off all, we assumed that all nuclear spins are of the 
same species. This is a commonly made approximation, 
where root-mean-square averages of the Zeeman and HI 
coupling constants can be used to mimic the situation of 
only one spin species being present. As long as the indi¬ 
vidual constants are not too different and if the sign of 
the nuclear magnetic moment is the same for all nuclear 
spins, we do not expect qualitative changes of our find¬ 
ings. Additionally, we have chosen all nuclear spins to be 
one-half. If one allows for larger nuclear spin quantum 
numbers I = 1/2, 3/2,..., the model can still be solved 
analytically. Since the thermal relaxation of the electron 
spin does only depend on its Zeeman energy, the tem¬ 
perature T\ is unchanged. In contrast to this, the spin 
flip temperature r 0 = p~ l AI(I + l)/3 increases by an 
/-dependent factor. Finally, we also implicitly assumed 
that all nuclei carry a spin. Yet, this is not the case for 
all materials. If I\ out of N nuclei carry a spin, the prob¬ 
ability to find the electron at the site of a nuclear spin 
is \(f>{fk)\ 2 = iV _1 = A' -1 • ( K/N ) = AT -1 • nj, where 
0 < nj < 1 is the abundance of spin carrying nuclei) 27 


Having convinced ourselves, that the results obtained 
within the simplified model should be reasonable for real 
systems, we finally want to give several examples, where 
we expect the spin flip to occur. The most severe con¬ 
straint is the negative sign of g *. As can be seen from 
Tab. [TJ this is realized, for instance, in III-V heterostruc¬ 
tures, where the electron experiences a strong spin-orbit 
interaction. Most promising among the considered ma¬ 
terials is GaAs, since its spin flip temperature is on the 
order of mK. The other materials having a smaller To suf¬ 
fer mostly from a large g* factor. Beside a small g* factor, 
potential materials would also benefit from a strong HI 
and from heavy nuclei with large g^ factors and large 
spin quantum numbers I. Among them, also systems 
with a negative gw such as CdTe can be considered, since 
this sign changes both the nuclear Zeeman coupling and 
the sign of the HI. Redefining the nuclear spin operator 
by I z -+ — Iz then yields the same results. 

Finally, we will briefly discuss the nature of the spin 
flip. If one leaves the equilibrium state while heating up 
the system (non-adiabatically), it will take some time for 
the system to reach the equilibrium at its new tempera¬ 
ture. Especially for crossing T 0 , the system is not only 
forced to flip the electron spin, but also approximately up 
to k > 1 nuclear spins. Thus, depending, for instance, 
on the microscopic details of the coupling of the spins to 
one or several baths, the time needed to equilibrate could 
be comparably long. 
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Appendix A: Diagonalization of the full Hamiltonian 

In this section, we show how the full Hamiltonian 
H a = H(a,v,a) in Eq. ^ can be diagonalized. 22 25 
Like the simplified Hamiltonian H 0 , it is invariant 
under an exchange of two nuclear spins. Moreover, 
a total flip of all spins (F : F\ms, rtiK, ■ ■ ■, mi) = 

| — ms, —itik, •■., —mi)) yields FH(a,u,a)F t = 
H(—a,—t/,a). Finally, the Hamiltonian commutes with 
the 2 -component of the total spin S z + ./- j 22 l 2! l Due 
to this property, it is convenient to represent the full 
Hamiltonian in the basis of product states between 
the electron spin and the total nuclear spin | ms) ® 

| J,Mj,{q i }) = | ms,J,Mj,{qi}), in which it has a sim¬ 
ple block-diagonal structure. The additional quantum 
numbers {qi} are related to the corresponding Clebsch- 
Gordon coefficient J —11—' l These quantum numbers q-; = 
qi({Ik}k=i) depend on the quantum numbers Ik of the 
original states | Ik,mk)- Since the Hamiltonian is degen¬ 
erate in these quantum numbers, we will use the short¬ 
hand notation J, Mj) for the nuclear spin state when¬ 
ever it is appropriate. Ordering these product states 
| ms, J, Mj) properly, this yields a block-diagonal rep¬ 
resentation of the Hamiltonian, where the blocks are of 
dimension 1 or 2. 

In this basis, all diagonal entries of the Hamiltonian 
are formed by the energies 

Ei(a, v, Mj, m s ) = {m s , J, Mj\H(a, v, a)\m s , J, Mj) 

= —ams — vMj + 2 msMj , (Al) 

where ms = 1/2,—1/2. The states 11/2, J, J) and | — 
1/2, J, — J) are already eigenstates of the Hamiltonian, 
and, hence the corresponding energies constitute the one- 
dimensional blocks. The off-diagonal parts of the two- 
dimensional blocks are given by 

F(J,Mj ) = (~,J,Mj\H(v,u,a)\±,J,Mj - 1) 

= i/J(J + 1) — Mj(Mj — 1) (A2) 

Thus, the eigenenergies and eigenstates of the two- 
dimensional blocks are obtained by diagonalizing the 2x2 












matrices 


H 2 (a, v,a,J,Mj) = 

E\{<j, v, M/, — |) aF(J,Mj) 
aF(J,Mj) Ei(a, v, Mj — 1, |) 


(A3) 


The eigenenergies of the 2x2 matrix in Eq. (A3) are, 
then, given by 



E 2 ,±(a, v,a,J,Mj) = 

^[Ei(a, v, — +£ , i(cr, u,Mj - 1, ^)± 

{ [Ei (u, - Ei (u, i/,Mj - 1, i)] 2 

+ 4 a 2 F(J,Mj) 2 }i (A 4 ) 


FIG. 4: (Color online) Dependence of the thermal expectation 
value of the electron spin (S Z ) T at t = 0 on the number 
of nuclear spin with (blue circles) and without (red squares) 
the flip-flop terms. The remaining parameters are p = 0.35, 
v = 2, and a = Kpu. 


Appendix B: Calculation of the partition function 

In this section, we want to calculate the partition 
function Z = Tr[e^' 3fl “] of the full Hamiltonian H a = 
H(<j, v, a) in Eq. (Jbj). Since the Hamiltonian is invariant 
under an exchange of two nuclear spins, the trace can be 
written as 


replacement ! — 11—11 is given by 

I m s ,'&M J )= H c fc 9i} l TO s. y -k, M j,{qi}) 

fe =0 {<Ji} 

(B4) 

where the coefficients obey the following relation 


Z = Y1 N Mj( m s,^Mj\e \m s , ^Mj), 

m s Mj=-% 

(Bl) 


where the nuclear spin state is given by 

IT \ I 1 11 1\ 

I ^Mj) = | • • ■ >~2> 2’“ " ’ 2 


(B2) 


--f — Mj ~2~ -\-ivij 


+Mj 


and = (k_ m )• In order to benefit from previ¬ 

ous resu itpmf 

which have been obtained for Mj > 0, 
we additionally use the effect of a total flip of all spins 
F\ms,^Mj) = I - rns^-Mj) and FH(a,u,a)F^ = 
H(—<j,—v,ol). Thus, the partition function can be ex¬ 
pressed by 

Z = 

K_ 

m s Mj= k mod 2 


N Mj(ms^Mj\e ljH{ a ’ v ' a) \ms^Mj)]- 

(B3) 


Mj — 


1 -K mod 2 


Finally, we replace the states | ms,'&M J ) by the 
states | ms,J,Mj). For spin one-half particles, this 


/ , \ u k 
{?,} 


(t - Mj)\(y + Mj)\{K — 2fc + 1) 
(K — k)\k\(K — k + 1) 


= d k 


(B5) 

The states 11/2, Mj, Mj, {%}} are already eigenstates of 
the Hamiltonian. Moreover, the state | — 1/2,0,0, {q,}) 
is also an eigenstate. All other remaining states can be 
expressed in terms of eigenstates of the Hamiltonian ob¬ 
tained from diagonalizing # 2 ( 0 ", u, a, J, Mj) in Eq. (A3). 
With this, we have calculated the partition function 
Z = Z(cr, v, a , K ), which allows us to find the thermal ex¬ 
pectation value of the electron spin according to Eq (11). 


In order to confirm that there is still a spin flip in the 
presence of the flip-flop terms, we calculated this thermal 
expectation for up to I\ = 60 nuclear spins. In Fig. [4j 
we show the electron spin at zero temperature. Clearly, 
the flip-flop terms reduce the polarization of the electron 
spin for small system sizes. For larger K , the thermal ex¬ 
pectation value (S z ) 0 , however, tends to the same value 
as without the flip-flop terms. In the next section, we 
analytically show, that (S z )0 = —1/2 is exactly reached 
for large system sizes. Moreover, we have analyzed the 
maximum value of the electron spin (S z ) 2 /n, which is ap¬ 
proximately found at r = 2 p 1 c.f. Fig. [ij As is clear 
from Fig. tal this quantity tends to {S z ) 2 / p = 1/2 for 
a large number of nuclear spins K. Note, that also for 
a = 0, this quantity is a function of K. Since the quan¬ 
tum fluctuations due to the flip-flop terms are already 
at zero temperature irrelevant for large system sizes, it 
is evident that they are also unimportant at higher tem¬ 
peratures. 
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FIG. 5: (Color online) Dependence of the thermal expectation 
value of the electron spin ( S Z ) T at r = 2 /p on the number of 
nuclear spin with (blue circles) and without (red squares) the 
flip-flop terms. At this temperature the electron reaches its 
maximum. The remaining parameters are p = 0.35, v = 2, 
and a = Kpv. 


Appendix C: Ground state of the Hamiltonian 

In this section, we want to analytically find the ground 
state of the full Hamiltonian H\. We choose er = Kpv , 
p < 1 and pv < 1, which are the necessary conditions for 
the spin flip. Moreover, we assume the Zeeman energy 
of a single nuclear spin to be v > 1, which can be always 
reached by a suitable choice of the external magnetic field 
in an experiment. In order to find the ground state en¬ 
ergy, we have first identified the minimum of the ener¬ 
gies E 2 -{<j , v, 1, J, Mj) for a fixed J with respect to Mj. 


From this set of minima, we have then found the smallest 
energy with respect to J. Finally, we compare this min¬ 
imum to the energies Ei(a, v, J, mg) resulting from the 
one-dimensional sub-spaces of the Hamiltonian. Follow¬ 
ing this scheme, we have found the ground state e nergy 
to be given by E 2 _(a, v, 1, K , K) defined in Eq. (A4|. 
The corresponding ground state is given by 


\E 2 - ( Kpv , v, a, J, Mj)) = 

V- (Kpv, v,a,J,Mj )| i, J, Mj - 1) 

+ (Kpv, v, a, J, Mj)\ - J,Mj), (Cl) 


where 


, K I< r-s 
r]_ 2 {Kpv,v,l, —, —) = 


2 2 ' t [(^)2 + l]i 


(C2) 


and 


i K K 

V-(Kpv,v,l, —, —) = 


2 2 [(^)2 + 1 ]i 


(C3) 


with r = 1 — K — v + Kpv , t = 2 K 2 and, s = (r 2 +1 2 ) 2 . 
In the limit of large K , {Kpv, v, 1,1^, K) = 1 an d 

rp! 2 (Kpv, v, 1, K , K ) = 0 and, hence, also for the full 
Hamiltonian Hi the ground state exhibits ( S z ) 0 = —1/2. 







